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Abstract Following up on hints of anisotropy in the cosmic microwave back- 
ground radiation (CMB) data, we investigate locally rotational symmetric 
(LRS) Bianchi type I spacetimes with non-minimally coupled scalar fields. 
To single out potentially more interesting solutions, we search for Noether 
symmetry in this system. We then specialize to the Brans-Dicke (BD) field in 
such a way that the Lagrangian becomes degenerate (nontrivial) and solve the 
equations for Noether symmetry and the potential that allows it. Then we find 
the exact solutions of the equations of motion in terms of three parameters and 
an arbitrary function. We illustrate with families of examples designed to be 
generalizations of the well-known power-expansion, exponential expansion and 
Big Rip models in the Friedmann-Robertson- Walker (FRW) framework. The 
solutions display surprising variation, a large subset of which features late-time 
acceleration as is usually ascribed to dark energy (phantom or quintensence) , 
and is consistent with observational data. 

Keywords Bianchi type I spacetime; Noether symmetry; Brans-Dicke theory. 



1 Introduction 

The recent decade has witnessed the observational evidence for the cosmic 
acceleration of the universe, which has become the central theme of modern 
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cosmology. This evidence is built up of observations of supernovae Type la 
(SNe la) [T], cosmic microwave background radiation (CMB) |2J, and large- 
scale structure of the Universe [3] . In Einstein's gravity, this acceleration can- 
not be explained by normal matter or fields on Friedmann-Robertson- Walker 
(FRW) metric background; therefore a mysterious cosmic fluid with negative 
pressure, the so called dark energy, is introduced. Presently, understanding 
the nature of the dark energy is one of the main problems in the research area 
of both theoretical physics and cosmology. A simple example of dark energy 
is the cosmological constant, equivalent to a fluid with the equation of state 
(EoS) parameter w = — 1, where w — p/p in which p is the pressure of dark 
energy, and p its energy density. However, the cosmological constant model is 
subject to the so-called fine-tuning and coincidence problems [I]. 

Dark energy, if it is a perfect fluid, must have EoS parameter w < —1/3. If 
— 1/3 > w > —1, the dark energy is called the quintessence [5]; if w < — 1, dark 
energy is dubbed a phantom fluid [5] . In spite of the fact that these models 
violate both the strong energy condition p + 3p > and the dominant energy 
condition p + p > 0, and therefore may be physically considered undesirable, 
the phantom energy is found to be compatible with current data from SNe la 
observations, CMB anisotropy and the Sloan Digital Sky Survey (SDSS) [TJIU 
01- 

Dark energy behavior can be exhibited by possible new fundamental fields 
acting on cosmological scales. The simplest such model is the single component 
scalar field. Other alternatives in the literature include the Cardassian expan- 
sion scenario [7], the tachyon [5], the quintom [5], the k-essence [TU| models 
and more. 

Alternatively, modifications of Einstein's gravity have been proposed to 
explain the cosmic acceleration of the universe; among them f(R) theories 
[llj . For example, quintessential behavior of the parameter w can be achieved 
in a geometrical way in higher order theories of gravity |12j . Scalar fields with 
various couplings and potentials can be put in by hand, or follow from the 
model naturally. 

One of the simplest modification of Einstein's gravity is the Brans-Dicke 
(BD) gravity theory [13j . a well known example of a scalar-tensor theory which 
represents the gravitational interaction using a scalar field in addition to the 
metric field. This theory is parametrized by one extra constant parameter, W. 
In the limit W — ► oo, the BD theory reduces to Einstein's gravity [TJ]. The 
conditions, where the dynamics of a self-interacting BD field can account for 
the accelerated expansion, have been considered in Ref. |16j , where it was con- 
cluded that accelerated expanding solutions can be obtained with a quadratic 
self-coupling of the BD field and a negative equation of state (EoS) parame- 
ter. Astrophysical data indicate that EoS parameter w lies in an interval of 
negative values roughly centered around — 1. 

The majority of popular cosmological models, including all the ones re- 
ferred to above, use the cosmological principle, that is, they assume that the 
universe is homogeneous and isotropic. On the other hand, there are hints in 
the CMB temperature anisotropy studies that suggest that the assumption of 
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statistical isotropy is broken on the largest angular scales, leading to some in- 
triguing anomalies |17j . To provide predictions for the CMB anisotropies, one 
may consider the homogeneous but anisotropic cosmologies known as Bianchi 
type spacetimes, which include the isotropic and homogeneous FRW models 
[15] . In this study we consider the simplest of these, the locally rotationally 
symmetric (LRS) Bianchi type I spacetime as an anisotropic background uni- 
verse model; note that this spacetime is a generalization of flat (k = 0) FRW 
metric. Our aim is to investigate the solutions of the field equations of scalar- 
tensor, in particular, BD theory of gravity for the LRS Bianchi type I spacetime 
using the Nocther symmetry approach. 

The Noether symmetry approach was introduced by De Ritis et al. [Ti?ll20] 
and Capozziello et al. [2T1I22] to find preferred solutions of the field equations 
and the dynamical conserved quantity. The Noether theorem states that if the 
Lie derivative of a given Lagrangian L dragging along a vector field X vanishes 



then X is a symmetry for the dynamics, and it generates a conserved current. 
Recently some exact solutions have been presented in the scalar tensor theo- 
ries following the Noether symmetry approach that allows the potential to be 
chosen dynamically, restricting the arbitrariness in a suitable way [T9ll20ll2Tl 



This paper is organized as follows. In the section [3J we present the field 
equations in scalar-tensor theory for the LRS Bianchi type I spacetime. In 
section [3l we search the Noether symmetry of the Lagrangian of scalar-tensor 
theory for the LRS Bianchi type I spacetime, and for the BD case, find new 
variables that include the cyclic one. In section @] we derive the EoS of the 
BD field. In section [SJ we give the solutions of the field equations by using 
new variables obtained in section [31 In section [BJ we choose some examples 
corresponding to the well-known solutions in the isotropic cases such as power 
law, exponentially expanding and phantom/Big Rip models. We also make a 
conformal transform to the so-called Einstein frame in Section [7] and discuss 
the problem in that context. Finally, in section [SI we conclude with a brief 
summary. 

2 The Lagrangian and the field equations 

The general form of the action that involves gravity non-minimally coupled 
with a scalar field is given by [2S], such that 



£*L = 0. 



(1) 



miiiiMiiisiissjiiziiiHiiiniissi . 




W{<P) 




(2) 



Here R is the Ricci scalar, F(<P) and W(<5) are generic functions that describe 
the coupling, t/(<?) is the potential for the scalar field <P, and <P a = <P- a stand 
for the components of the gradient of the scalar field. Note that we use Planck 
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units. For F{<&) = 1/2 and W(<5) = #/2, the action reduces to the form of 
Einstcin-Hilbert action minimally coupled with a scalar field. The choice of 
F(<P) and W(<P) give us other gravity theories such as the BD theory, for 
which F((P) = <P and W(<P) =constant. 

Variation of the general form of the action with respect to metric tensor 
yields the field equations 



F($)G ab = T* 



where G a b = Rab — hR9ab is the Einstein tensor, 



rp<P 

1 ab 



W 



$n$b 



W 



9ab$c$ c ~ ^9abU{<P) - g ab aF(<P) + F(<P). ab 



(4) 



is the energy-momentum tensor of scalar field, and □ is the d'Alembert op- 
erator. It is clear that Tf b includes the contributions from the non-minimal 
coupling and the scalar field parts of the action varying with respect to metric. 
In Ref . [31| , it has been discussed that there are three possible and inequivalent 
ways of writing the field equations, corresponding to the ambiguity in the def- 
inition of the energy-momentum tensor of non-minimally coupled scalar field. 
The variation with respect to <P gives rise to the generalized Klein-Gordon 
equation governing the dynamics of the scalar field 



W 

2 — D$ + RF'($) 



w 



w 

<£2 



<P C <P C - U'{<P) = 0, 



(5) 



where the prime indicates the derivative with respect to <P. Note that this 
equation is equivalent to the contracted Bianchi identity. It follows from G a b b — 
that using Eq.© together with Eq. ©, yields 



1 

F' 



-X 



&(*) 



= 0<=>T 



b(t) 



a ;b 



p a 



(6) 



From ©, we compute T a b ^ to get 



b(t>) 



W\ W 1 



F.bRn 



(7) 



Using Eq.©, Eq.© and F. b = F'($)$ b in Eq.©, we have 

'W 



w 

2—D$ + RF'($) 



(p <p2 



<P C <P C - U'($) 



= 0. 



(8) 



This is obviously the generalized Klein-Gordon equation for <P a ^ 0. 
The line element of the LRS Bianchi type I spacctimc has the form 



ds 2 = -dt 2 + A 2 dx 2 + B 2 (dy 2 + dz 2 ) , 



(9) 
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describing an anisotropic universe with equal expansion rate in two of the 
three dimensions. The Ricci scalar of this spacetime is 



R 



B 
'B 



B*_ 
B 2 



AB 
'AB 



(10) 



where the dot represents differentiation with respect to t. For the metric 
the field and generalized Klein-Gordon equations can be obtained from Eqs.© 
and ([5|) respectively 



B 2 AB 

h 2 

B 2 AB 

B B 2 

2 1 h 

B B 2 
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2F 
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<P 2 (W W 
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U' 
2F 1 
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= 0, 
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2F 
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F 7 ® 









0. 



(11) 

(12) 



(13) 



(14) 



where F' ^ 0. The Lagrangian density of the LRS Bianchi type I spacetime is 



L = -2FAB 2 - AFBAB - 2F'B 2 A$ - 4F'ABB& 
W<P 2 



-AB 



(15) 



Using such a Lagrangian, one may obtain the Eulcr-Lagrange equations as 
given in (|T2"1) - (|14I) . The energy function, El, associated with the Lagrangian 
(T13]) is found as 



dL • dL ■ 8L- 

E L = —?A H rB H r<P 

dA dB d<P 



B 2 



AB 
'AB 



F' 








— <P 2 + U 






< 




F 






2F 





(16) 



Therefore, it is obvious that the (0,0)-field equation given by ([TTj) is equivalent 
to E L = 0. 
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3 Noether symmetry approach 

In this section we seek the Noether symmetry of the the Lagrangian (|T5|) . The 
configuration space of this Lagrangian is Q = (A, B, <?), whose tangent space 
is TQ = (A, B,<P, A, B,<P). The existence of Noether symmetry implies the 
existence of a vector field X such that 

v 9 n 9 8.8-8.8 
X = q— +/?— +7— + +7^ (17) 

OA ' 8B 'd$ dA dB 8<P 

where a, and 7 arc depend on A, B and <P. Hence the Noether equation given 
by © yields the following set of equations 

80 F' di 

2 M + B T^">- < 18 » 

a „ da ,80 ,F'fj^8j\ n .„„. 
- + B— + A-- + A— [ - + B—-) = 0, (19 
2 8B dB F \2 dB J ' v ' 

2 4 + f{^ + B 8A +B 8^ + 2A £ 
F" W di 

+ ^-F* AB 8A=°> ™ 
da A 80 F' ( A B da 80 <9 7 

d^ + Bd^ + T\ a + W + 1 2dB +A dB +A d$ 

^ WAB^ = 

F ' 2F<P dB v ' 

(Ba + 2A0)U + AB jU' = 0. (24) 

Leaving these equations as reference for future work, we now specialize to the 
Brans-Dicke case, where F = <P and W is constant. We also require the Hessian 





fW 











determinant, D 



to vanish in order to get nontrivial solutions. 



This condition reads for our case as 

1 (\ A W 

D = — — (3<Z>F' 2 + 2WF) = 0, (25) 

which also determines a W value. So we will be working with the coupling 
functions 

F = $, W = ~. (26) 
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Thus the Lagrangian (|15|) becomes degenerate, and the BD action has the 
form 



<PR + — <P C <P C - U($) 



(27) 



which can be easily related the conformal relativity by defining new scalar 
field cp as & = </? 2 /12, transforming the action into 



d Xy/~g 



<p 2 R + <p c <p c - U(ip) 



(28) 



This action is of course conformally invariant, since the application of the 
conformal transformation formulas together with the appropriate integration 
of the boundary term gives the same form of this action (see Refs. [Tl] and 
[27]). 

Dabrowski et al. |27j have shown that the anisotropic non-zero spatial cur- 
vature models of Bianchi types I, III and Kantowski-Sachs type are admissible 
in W = —3/2 BD theory. They have given solutions of BD field equations 
for these spacetimes without the potential U(<P). In a previous work [28], we 
studied the case F = -^((p — <I>o) 2 i W = ^ for the Bianchi types I and III, and 
Kantowski-Sachs spacetimes, where e is a parameter depending the signature 
of metric, but did not study the qualitative (acceleration etc.) behaviour of 
the universes in the solutions found. 

In our present study we find the solutions of BD field equations with po- 
tential U(<P) for the LRS Bianchi type I spacetime, although we derive the 
equations for a general function W($). In this context, the solutions of the 
above set of differential equations (|18l) - ([2~i)) for a, (3, 7 and potential U($) are 
obtained as 



\<P 2 



I ( A 2 $ 3 ^ 



A 2 B<P 1/2 



(29) 



where A is a constant. The potential found indicates that for the Noether 
symmetry to be present, the scalar field must be a (massive) free field. The 
vector field X generating the Noether symmetry and determining the dynamics 
of LRS Bianchi type I metric is given by 



X 
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4 Equation of state 

Using F(<P) — <P and W = —3/2 in Eq. Q we can derive the energy density 
of the scalar field 

^_^ + m_ 3[m (31) 

and the directional pressure due to the scalar field 
3 2 U($) 

pf = —_- -Li +* + (3ff i = x,y,z (32) 

for the x, y and z directions. Here Hi represents the directional Hubble pa- 
rameters in the directions of x, y and z respectively, and may be defined as 



H x — — , H y = H z - —. (33) 



H = - = - (H x + 2H y ) , (34) 



— , H<t, — H z — — . 
A B 

The mean Hubble parameter for LRS Bianchi type I metric is given by 

a _ 1 

a 3 

where we define a = (AB 2 ) 1/3 as the average scale factor of the Universe. Also 
of interest is the volumetric deceleration parameter (for the isotropic case, 
simply the " deceleration parameter" ) , denned by 

q = -^ = - l -m (35) 

so that it is dimensionless. 

The pressure is a vectorial quantity, as would be expected for anisotropic 
expansion, and thus the EoS parameter of the scalar field for the LRS Bianchi 
type I spacetime may be determined separately on each spatial axis: wf(t) — 
pf/p*, i.e. 

and the average EoS parameter of scalar field for the LRS Bianchi type I 
spacetime may be defined as 

1 2£i + EM _<z>_ 2H& 

w* = ~* = ~ K + 2w y ) = 4g> .T 2 ,u,, (37) 

where the is the isotropic pressure |32j . 

In the anisotropic case, w does not directly determine the sign of q as it 
does in the well-known (spatially) flat FRW cosmologies. To understand this, 
consider the Raychaudhuri Equation |33j 

^ + l^ 2 + ~ + \ (P + 3P) ~ A = ° (38) 
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where 6 is the expansion scalar. a^ v is the shear tensor, uj^ is the vortic- 
ity tensor (descriptions/definitions also given in [33]) an d we take A = 0. If 
comoving particles in a metric move geodesically, the expansion scalar also 
gives the expansion rate of the space itself, when we apply this equation to a 
swarm of such particles. This is the case for the LRS Bianchi I metric, where 9 
turns out to be equal to 3H, therefore 8 + ^9 2 is —3H 2 q. The vorticity tensor 
vanishes for this case, and 



2B{t) 2 
3A(t) 2 



d (A(t)\ 
dt \B{t) J 



(39) 



Therefore in the isotropic case a^a^ also vanishes and since in that case 
p cx H 2 , we find that q is proportional to 1 + 3w. This correspondence is 
obviously broken when there is anisotropy, and therefore w is not that useful 
a parameter. 



5 The solutions from new coordinates and Lagrangian 

The solutions of dynamical Eas. ([TTj) - ([l"4"|) are not easy to evaluate in the present 
form. In order to simplify these equations we can search for the cyclic vari- 
able (s). In case of Noether symmetry, we should introduce new variables in- 
stead of old variables, i.e. a point transformation {A, B,<P\ — > {/i, v, u} in 
which it is assumed that fi is the cyclic coordinate. A general discussion of 
this procedure has been given in Ref . |34j . The new variables {/i, v, u} satisfy 
the following equations 

iy^dn = 1, iy^dv = 0, iy^du = (40) 



where ix is the interior product operator of X. Using Ea. ((29|) . a solution of 
the above Eqs.gO]) yield 



f i = A 2 B$ 3/2 , v = B$ 1/2 , u = A<P. (41) 
The inverse transformation of these variables are 

A = —, B= { -^, (42) 



from which it follows that 



i{t) - (43) 



1 /3 

for the average scale factor of the Universe defined by a — (AB 2 ) . Con- 
sidering the transformation of variables, the coupling functions (|2^|) and the 
potential U('P) given in ([2^1) . the Lagrangian (TT5j) becomes 

L = -lu^iiv - Xuv 2 (44) 
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which does not depend on /i (i.e. ^ = 0), as desired. This independence is 
preserved if we make further transformations in the u — v plane; and also is 
inherent in the procedure described in Ref.[34]. We used this freedom to get 
as simple a Lagrangian as possible. 

The Euler-Lagrange equations relative to this Lagrangian are 



v — Iqu = 
A 

A" - 2- 
fill 



2 2 
V U 



0, 



= 0, 



(45) 
(46) 

(47) 



where Iq is a constant of motion associated with the coordinate fi. Here ix@L = 
ao = —21q and Ol is the Cartan one-form, and Ea. (|46[) is equivalent to the 
vanishing of the energy functional, as noted before. 
By substituting u from Ea. (j45j) into Ea. (|46|) we get 



fi{t) = l x v 3 + l 2 



(48) 



where l\ = A/(6/q), lo ^ 0, and I2 is an integration constant. The remaining 
equation (|47|) is identically satisfied now. Inserting (|45l) and (|48|) into (|42|) . 
we obtain the metric functions and scalar field in terms of the arbitrarily 
specifiable v (henceforth to be called the seed junction), and u = v/l® as 
follows: 



A(t) = h 



B(t) 
#(*) = 



lo 

m 

Ht)Ht) 2 

l 2 (lMt) 3 +l2)' 



Thus, using P^|) and (|50p . the LRS Bianchi type I metric becomes 



-dt 2 + ll (hv{tf + l 2 ) 



m 2 



h + 



dx" + dy + dz* 



Then, the average scale factor of the Universe takes the form 

ait) = ^ {hvitf + hf 3 . 



(49) 
(50) 
(51) 

(52) 

(53) 



Also the case of isotropy can be seen, by requiring A oc £?, to correspond to 
h = 0. 

In the general anisotropic case, the directional Hubble parameter in x di- 
rection defined in Eq. ([3"3"l) is found as 





V 


u 


V 


V 


3Zi^ 3 


fj, 


V 


u 


-- - 

V 


V 


his 3 +l 2 



- 1 



(54) 



Noether-symmetric LRS Bianchi type I universes in scalar-tensor BD theory 



11 



the other directional Hubble parameters in the directions y and z are similarly 
found as 



u v v 
u v 2v 



_hv 3 + h 

and the mean Hubble parameter defined in Eq. (|34p is found as 

2hv 2 v 



+ 1 



3 [i u 



(55) 



(56) 



At this point, we would like to emphasize that particular solutions can be 
found, that is, the scale factors A(t),B(t), the potential <P(t) and the average 
scale factor a(t) can be obtained by specifying the function v(t). In the next 
section, we will find such examples. 



6 Examples 

To illustrate our family of solutions for the dynamics of LRS Bianchi type 
I cosmologies containing a scalar field, we choose examples corresponding to 
the well-known solutions in the isotropic case, that is, power-law models, i.e. 
anisotropic t 71 , exponentially expanding models, i.e. a; so t r opic oc e x *, and phan- 
tom/Big Rip models, i.e. ai so t r opic oc (t c — t)~ CT , where 77, \ and a are positive 
constants. We can find the form of each v(t) by solving from Ea. (|53|) after 
setting l 2 = 0. 

The most general form we find for the power-law models is v (t) = C2 _^ v -i , 
but for simplicity we use v{t) = a\t n . For the exponential case, the most 
general form is v{£) = g4 _|_ xt , but we use u(t) = a2e kt . Similarly, for the 
phantom/Big Rip models, the most general form is v(t) = CB _^tty+i > but 
we use v(t) = , f °^ t \m (n, k or m are not necessarily integer). 



6.1 Power-law seed 



When we use v(t) = a\t n , where a\ and n are nonzero constants, the scale 
factors read 

A(t) = ^-^(halt 3 * + l 2 ), B(t) = -^-t 1 -^ 2 [ha\t Zn + l 2 ] 1/2 , (57) 
the average scale factor is found from Eq. ([5"3"|) as 

a(t) = -^- [halt^ + hfh 1 -*, (58) 



the scalar field from Eg. ([51]) as 
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and the mean Hubble parameter from Ea. (|34j) as 

2ha\nt^ 1-n 
H - l ia \t^ + l 2 + —- {60) 

Inspection of the scale factors A(t), B(t) and a(t) shows that if l\ or I2 
is zero, they will reduce to single powers of t. But also in the general case 
where both l\ and I2 are nonzero, the scale factors will approach one power of 
t near zero and another near infinities. These powers depends on the sign of 
n: For negative n, near t = we have A(t), B(t), a(t) oc t 1+n , near infinities 
we have A(t) oc t 1 ~ 2n , B(t) oc i 1- ™/ 2 and a(t) oc i 1_n ; for positive n, the 
behaviors near zero and infinities are interchanged. Therefore, we need to 
treat the cases of vanishing l\ or I2 separately from the general case. These 
cases will be split into subcases, for example, according to the value of n, etc. 
Each line in the Tables 11121 and [3] shows the timeline of an eventual subcase 
from t — > —00 to t — > +00, which usually describes more than one universe, 
since the timeline may be interrupted by singularities, shown by expressions in 
square brackets in the lines of the tables. Since solutions cannot be taken valid 
through singularities; each interval between singularities should be considered 
an independent solution/universe. The universes not extending to the right or 
left edge of a line have finite lifetime. 

Obvious candidates for singularities are t = and t — > ±00, in light of 
the above discussion. Moreover, the (Ziafi 3 ™ + I2) terms in the scale factors 
can vanish at some finite t value t\, or diverge at t = 0, if both l\ and I2 are 
nonzero. The meaning of a singularity not only depends on its mathematical 
behavior, but also on its relative time-relation to the observers: If observers see 
(or calculate) vanishing comoving volume in their finite past or finite future, 
they will call that singularity a Big Bang (BB) or a Big Crunch (BC); This is 
shown as "BC[0]BB" in the tables. Similarly, observers calculating diverging 
of the comoving volume in their finite future will call that singularity a Big 
Rip (BR), and the case of diverging comoving volume in the finite past we will 
call it an Inverse Big Rip (iBR). 

The square brackets should contain an ordered triple, where the first entry 
shows the behavior of A(t), the second entry the behavior of B(t) and the third, 
the behavior of a(t). The behavior is indicated by the symbols 0,C or 00 to 
denote that the relevant scale factor vanishes/goes to a finite number/diverges 
at that point. To make the tables more compact, an ordered pair (describing 
the behavior of A(t) and B(t); whenever these make the behavior of a(t) 
obvious) or a single entry (indicating that the behavior of all three scale factors 
is the same) may be used. Finally, one should also note that negative t can 
only be considered if n is a rational number with an odd denominator, not for 
general n. 

We discuss the three above-mentioned cases in order of increasing com- 
plexity: First, we discuss the isotropic I2 — case, then the simple anisotropic 
l\ = case, and finally the general anisotropic case where both l\ and I2 are 
nonzero. 
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6.1.1 The I2 = (isotropic) case. 

As stated after Eq. ([55| , in this case the LRS Bianchi type I spacetime becomes 
isotropic, i.e. reduces to the well-known spatially fiat FRW metric. The scale 
factors become proportional to each other, and 

a{t) = a Q t 1+n (61) 

if/ 3 

where clq = fli- 5 -^ — ■ This means that a(t) can vanish or diverge at zero or 
infinity, and the critical n values are -1 and 0, leading to the subcases shown 
in Table [T] The cases are 



t 

n 


—00 









+00 


n < -1 


[0]Z 




BR[oo]iBR 




Z[0] 


n= -1 


(C) 


— static — 


(C) 


— static — 


(C) 


-1 < n < 


[oo]I 


"decelerating" 


BC[0]BB 


decelerating 


I[oo] 


n ->• 


[oo]I 


linear 


BC[0]BB 


linear 


I[oo] 


n > 


[oo]I 


" accelerating" 


BC[0]BB 


accelerating 


I[ooj 



Table 1 Universes derived for I2 = (isotropic universes) where v(t) is a power of t. Each 
line shows a timeline from —00 to 00, but may represent multiple universes delimited by 
singularities, shown by square brackets. The expressions in the brackets refer to the scale 
factor a(t). Also, I: infinite scale factor at infinite future or past, BR: Big Rip, i: Inverse, 
BC: Big Crunch, BB: Big Bang, Z: vanishing scale factor at infinite future or past [Also see 
text and the note about negative time on pagc ll2l . 



• n < — 1: Negative time£] describe a universe that expands from a van- 
ishingly small scale factor in the infinite past, and then the scale factor 
diverges within a finite time: one might call this a Zero-Big Rip (Z-BR) 
universe. Positive times describe another universe, whose behavior is the 
time-revcrse of the first, an iBR-Z universe. 

• n = — 1: This universe, being static and flat, is equivalent to Minkowski 
space. 

• — 1 < n < 0: For positive times, this case represents decelerating (spa- 
tially) flat FRW universes, including the radiation-dominated (?i = —1/2) 
and matter-dominated (n = —1/3) expansions, in particular. For negative 
times 1 , however, the universes contract from infinite scale factor in the in- 
finite past to a Big Crunch at t = 0; so the universes represented on the 
last three lines of Tabic Q] may be called Infinity-Big Crunch (I-BC) and 
Big Bang-infinity (BB-I) universes. The I-BC universes are decelerating, 
too: The negative a(t) [for positive a{t)] becomes more negative with time. 

• n — > 0: This case can be analyzed with the help of an infinite scale trans- 
formation [note that Eqs. (|45|) - (|47)) arc satisfied for constant u(t)], and 
describes universes linearly expanding (positive time) or contracting (neg- 
ative time 1 ) with increasing time. 



1 subject to the condition on n mentioned on page 1121 
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• n > 0: Positive times describe a universe that starts with a Big Bang and 
expands forever with increasing speed; negative times 1 describe a universe 
with time-reversed behavior. 



As is well-known (and mentioned above) , in FRW cosmology the behavior 
of the universe and the properties of the (possibly efective) fluid contained in 
the universe are related. In fact in our I2 = case the Hubble parameter, the 
deceleration and EoS parameters and the scalar field simplify to 



H 



n + 1 



q = 



*(*) 



n + l' 
3n+ f 

"3(n+ 1)' 

6n 2 



A 



-r 



(62) 
(63) 

(64) 
(65) 



and we can confirm that q is proportional to 1 + 3U;; in particular, for n = —1/2 
Ea. (|64)) gives w = 1/3 and for n = —1/3, it gives w = 0. For n — >• 0, the 
deceleration parameter vanishes. We know that empty universes can have this 
property, and we see that the scalar field vanishes in this case. For n > 0, 
the effective fluid is equivalent to quintessence/quiessence (— 1 < w < —1/3), 
approaching cosmological constant (w = — 1) as n — > 00. Only for n < — 1 do 
we have phantom behavior (w < —1), and this is the case where a Big Rip 
appears (for negative times 1 , at least). 

The n = — 1 case is particularly interesting. For this case, q and w seem to 
diverge, while in Table Q] this subcase does not seem to be more problematic 
than others. But, as n approaches —1, the scale factor a(t) approaches a con- 
stant function, therefore a(t) and d(t) approach zero. Hence the deceleration 
parameter, which has division by a 2 in its definition, is not a good parameter 
to use in this limit. Also, since this universe corresponds to Minkowski space, 
the stress-energy-momcntum tensor should vanish. Let us calculate the energy 
density and isotropic pressure, using Eq. (|65p and the potential expression in 
Eq.jinj), in expressions ([21]) and (pt2"j) : 



p* = Hn 2 (n+l) 2 i- 4 , 



(66) 
(67) 



These expressions verify that n — corresponds to empty universes. It can 
also be seen that both the density and pressure vanish for n = —1, although 
the field does not vanish! Of course, this also makes w meaningless in this 
subcase. 
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6.1.2 The h=0 case. 

Since l\ = A/(6Zq), in this case the potential f/(#) = A^ 2 vanishes, therefore 
this case is equivalent considering the scalar field to be massless. Now A{t) oc 
t 1 - 2 ", B(t) oc t 1 - 71 / 2 and a(t) oc t 1 "™, so that the critical n values are 1/2, 1 
and 2; and we are led to Table [2] listing the subcases. 



\- t 

n 


— oo 





+oo 


n < 1/2 


[oo, oo, oo]I 


BC[0, 0, 0]BB 


I[oo, oo, oo] 


n = 1/2 


[C, oo,oo]2I 


2BC[C, 0, 0]2BB 


2I[C,oo,oo] 


1/2 < n < 1 


[0, oo, oo]IP a 


cBC[oo, 0, 0]cBB 


IPa[0, OO, OO] 


n = 1 


[0, oo, C]IP b 


BD[oo, 0, C]iBD 


IP b [0, oo, C] 


1 < ra< 2 


[0, oo, 0]IP c 


cBR[oo, 0, oo]ciBR 


IP c [0, oo, 0] 


n = 2 


[0, C, 0]1Z 


lBR[oo, C, oo]liBR 


1Z[0, C, 0] 


n > 2 


[0, 0, 0]Z 


BR[oo, oo, oo]iBR 


Z[0, 0, 0] 



Table 2 Universes derived for l\ = 0, where v(t) is a power of t. Each line shows a timeline 
from — oo to oo, but may represent multiple universes delimited by singularities, shown by 
square brackets. The expressions in the brackets refer to A(t), B(t) and a(t). Also, I: infinite 
scale factors at infinite future or past, BR: Big Rip, i: Inverse, BC: Big Crunch, BB: Big 
Bang, 1 or 2: only one or two of the spatial dimensions, c: cigar-type, BD: Big Draw, IP a : 
"Infinite pancake" of type a, etc., Z: vanishing scale factors at infinite future or past [Also 
see text and the note about negative time on page !12l . 



Let us note that 

• For n = 1/2, we get universes evolving in the y and z dimensions only, 
either from infinite scale factors to a BC, or from a BB to infinite scale 
factors, in infinite time, possibly to be called 2I-2BC and 2BB-2I for two- 
dimensional 

• For 1/2 < n < 1, A(t) diverges at t = 0, while the other scale factors 
vanish. This singularity is extremely anisotropic: As the universe contracts 
towards a BC in the y and z directions, it actually expands infinitely in the 
x direction! But this expansion is not enough to prevent the vanishing of 
a comoving volume, and in this sense this event is a BB/BC. Singularities 
where one dimension diverges while two shrink to zero are called "cigar- 
type", so we might call this event a "cBC/cBB". 

Again for 1/2 < n < 1, A(t) vanishes as t — > ±oo, while the other scale 
factors do not. This type of singularity is also extremely anisotropic, and 
is called a "pancake" singularity. In this particular subcase, the other scale 
factors diverge, so the singularity is called an infinite pancake. This type of 
singularity appears in the next two subcases too, but the behavior of the 
average scale factor a(t) is diferent in each subcase. So we classified the 
"infinite pancakes" accordingly. (In [35], a classification is made according 
to the behavior of the three separate scale factors of a Bianchi I spacctimc.) 

• For n — 1, A(t) diverges, a(t) goes to a constant, B(t) vanishes at t = 0. 
Although this t — event is also a cigar-type singularity, since the comoving 
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volume neither vanishes nor diverges, it is not a BB/BC or BR/iBR. We 
suggest the names Big Draw (BDjj and Inverse Big Draw (iBD). 
• For n — 2, the universe evolves in the x dimension only, either from van- 
ishing scale factor to a BR, or from an iBR to vanishing scale factor, in 
infinite time. For n > 2, this occurs in all three dimensions. 

For li = 0, the mean Hubble parameter, the deceleration parameter, the 
EoS parameter and the scalar field simplify respectively to 

H= 1 —^, (68) 



t 
n 

— 5n - 



3(n-2)' 



(69) 
(70) 



*(*) = jn-t 3n ~ 2 - (71) 

There seem to be problems for n = 1 and n = 2, while in Table [5] these cases 
do not seem to be more problematic than others. The apparent divergence of 
q as n — 1 is understood as for the n = — 1 line of Table [T] a(t) becomes 
constant. As discussed at the end of Section |4l the w parameter is not very 
useful when there is anisotropjo, but one might still ask why it diverges for 
n = 2. Let us again calculate the energy density and isotropic pressure of 
the effective fluid for our subcase, recalling that the potential vanishes. Then, 
using Ea.([7Ijl. we get 



= J^(n-2)(3n-2)t 3 "-*, (72) 



f = -^-n 2 (2-5n)(3n-2)i 3 "- 4 . (73) 



explaining the divergence of w by the vanishing of for n = 2. We also see 
that the density and pressure both vanish (other than for the trivial, isotropic 
n = case) for n = 2/3, due to the constancy of the field and vanishing of the 
potential. 



6.1.3 The case with both l\ and I2 nonzero. 



For this subcase, the limiting behaviors of the scale factors mentioned on 
page IT21 show that they diverge near infinities for all n. Moreover, at late time, 
the universe expands with acceleration, again for all n, as can be seen from the 
asymptotic behavior mentioned after Eq. ((6"0"]) . But the behavior of the scale 
factors near zero determines the critical values of n as -I, 1/2, 1 and 2, leading 
to the subcases shown in Tabled] Outstanding features are 

2 since the volume of a metal drawn to produce a wire does not change. 

3 As a specific example, consider n = 7/5, which gives w = 25/9 according to H70H . which 
should give deceleration, but gives q = —7/2, that is, acceleration according to d 69 1) - 
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n, ti 


f. 


— oo 


-1*1 1 





1*1 1 


+oo 


n < — 1 


No ti 


[oo]I 


— B— 


BR[oo]iBR 


— B— 


I[oo] 




ti < 


[oojl 


BQOJBB 


BR[oojiBR 


— B— 


I[ooJ 




ti > 


|oo|I 


— B— 


BR[oo]iBR 


BQOJBB 


I[ooJ 


n = -1 


No ti 


[oo]I 




(C) 


— B— * 


I[oo] 




3*1 


|oo|I 




(C) 


BQOJ'BB 


I[oo] 


-1 < n< 


No ti 


[oo]I 




BQOJBB 




I[oo] 




*1 < 


[oojl 


BQOJBB 


BQOJBB 




I[ooJ 




ii > 


[oojl 




BQOJBB 


BQOJBB 


I[ooJ 


71 -> 


[oo]I 


linear contraction 


BQOJBB 


linear expansion 


I[oo] 


< n < 1/2 


No t\ 


[co]I 




BQOJBB 




I[oo] 




h < o 


[oojl 


BQOJBB 


BQOJBB 




I[ooJ 




ti > 


[oojl 




BQOJBB 


BQOJBB 


I[ooJ 


n = 1/2 


No ti 


[co]I 




2BQQ0J2BB 




I[oo] 




ti < 


[oojl 


BQOJBB 


2BQQ0J2BB 




I[ooJ 




h > 


[oojl 




2BQQ0J2BB 


BQOJBB 


I[ooJ 


1/2 < n < 1 


No ti 


[co]I 




cBC[oo,0,0]cBB 




I[oo] 




ti < 


[oojl 


BQOJBB 


cBQoo,0,0JcBB 




I[ooJ 




il > 


[oojl 




cBQoo,0,0JcBB 


BQOJBB 


I[ooJ 


n = 1 


No ti 


[co]I 




BD[oo,0,C]iBD 




I[oo] 




il < 


[oojl 


BQOJBB 


BD[oo,0,CJiBD 




I[ooJ 




h > 


[oojl 




BD[oo,0,CJiBD 


BQOJBB 


I[ooJ 


1 < n < 2 


No ti 


[oo]I 


— B— 


cBR[oo,0,oo]ciBR 


— B— 


I[oo] 




il < 


[oojl 


BQOJBB 


cBR]oo,0,ooJciBR 


— B— 


I[ooJ 




il > 


[oojl 


— B— 


cBR|oo,0,ooJciBR 


BQOJBB 


I[ooJ 


n = 2 


No ti 


[co]I 


— B— 


lBR[oo,C,oo]liBR 


— B— 


I[oo] 




il < 


[oojl 


BQOJBB 


lBR[oo,QooJliBR 


— B— 


I[coJ 




il > 


[oojl 


— B— 


lBR[oo,c,ooJliBR 


BQOJBB 


I[ooJ 


n > 2 


No ti 


[oo]I 


— B— 


BR[oo]iBR 


— B— 


I[oo] 




ti < 


[oojl 


BQOJBB 


BR[ooJiBR 


— B— 


I[ooJ 




h > 


[oojl 


— B— 


BRJooJiBR 


BQOJBB 


I[ooJ 



Table 3 Universes derived for l\ ^ and I2 7^ 0, where u(t) is a power of t. Each line shows 
a timeline from —00 to 00, but may represent multiple universes delimited by singularities, 
shown by square brackets. A single expression in the brackets refers to all scale factors, two 
expressions refer to A(t) and Bit), three expressions refer to A(t), B(t) and a(t). Also, I: 
infinite scale factors at infinite future or past, B: Bounce, BR: Big Rip, i: Inverse, BC: Big 
Crunch, BB: Big Bang, 1 or 2: only one or two of the spatial dimensions, c: cigar-type, BD: 
Big Draw [See text, and also note about negative time on page [12] *Sign of bounce time or 
t\ may be positive or negative] . 



• The comoving volume diverges at t = for n < — 1 and for n > 1. Since it 
also diverges as t — > ±00, this means that it must go through a minimum 
between two divergences — a Bounce (B). If t\ does not exist, there will 
be a Bounce each in the t < universe 1 and the t > universe; whereas 
if t\ does exist, it will essentially bring one of the minima down to zero, 
inserting a BC/BB, effectively splitting one of the universes into two. Hence 
the first line of Table |3] has an Infinity-Bouncc-Big Rip (I-B-BR) universe 
and an iBR-B-I universe; the second line an I-BC universe, a finite-lifetime 
BB-BC universe and an iBR-B-I universe; and so on. 
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• The third line features an iBR-BC universe, i.e. a finite-lifetime universe 
which starts with infinitely large scale factors that immediately contract 
to finite values, and keep on contracting to a BC. 

• For n = — 1, the scale factors become constant at t = 0. Therefore the sub- 
case without t\ features the only universe of the table without beginning 1 
or end, a bouncing universe. 

• The n — > case is the same as the corresponding case of Table [TJ In fact, 
it is isotropic. 

• For all cases with n > 0, the singularity at t = is qualitatively the same 
as in the corresponding n value Table [5J 

Interestingly, a wide range of cosmological possibilities, including late-time 
acceleration, can be reproduced by using a power of t as the seed function. This 
even includes the isotropic case, where the time-dependence of the source scalar 
field does not change with n: <P(t) oc t~ 2 . This surprising richness, including 
the possibility of vanishing energy-momentum tensor while the field is nonzero, 
seems to be due to the coupling coefficient selected by the Hessian determinant 
condition, and the potential selected by the Noether symmetry approach. 



6.2 Exponential seed 



When we use v(i) = a\e kt where a 2 and k are non-zero constants, in Eqs. (|49[) - 
(fSTj) . (|5"5)) . ([5B1 and (|35[) : the scale factors, the scalar field, the mean Hubble 
parameter and the deceleration parameter become 



A(t) = 
B{t) = 
a(t) = 

*(*) = 
H = k 



lo 



k^fai 
ka 2 



k 2 a\ 



(h4e 3kt + h) e- 2kt 

3„3fct 
2 e 

ikt 

3ki 



(ha 



(l\a 



3„3kt 



+ h) 1/2 e- kt / 2 , 
h) e , 



l 2 ha 3 2 e 3kt + W 
ha 3 2 e 3kt - 1 2 



ha 3 e 3kt + h 



2„6„6fet 



Uxha\e 3kt + 1\ 



(ha 3 e 3kt - l 2 y 



(74) 
(75) 
(76) 
(77) 
(78) 
(79) 



Unlike in Scction lBTTl these scale factors are never singular at t = 0, however 
the (Ziaje 3 *^ + 1 2 ) terms in the scale factors can vanish at some finite t value 
t 2 , if l\ and l 2 are both nonzero. Since they do not diverge at finite time, there 
is no Big Rip. 

In that case, the scale factors diverge exponentially at both infinities, re- 
gardless of the signs of l±, l 2 and k. Therefore, if t 2 does not exist, the universe 
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is a bouncing universe, infinite in both time directions. This solution is quali- 
tatively similar to the one in the fourth line of Table |3J If £2 does exist, then 
the solution represents two universes; one collapsing from infinite scale factors 
to a Big Crunch, one expanding from a Big Bang to infinite scale factors. This 
solution is qualitatively similar to the one in the fifth line of Table [31 For all 
solutions (with nonzero l\ and I2), at late time, H — > \k\ and q — > —I, con- 
firming the de Sitter- like asymptotic behavior of the universe(s). For positive 
k, the scalar field ^ approaches a constant. 

If either one of l\ and I2 is zero, the universe either expands or contracts 
exponentially, depending on the sign of k and on which constant vanishes. In 
these cases, q — —I at all times, and for I2 — (isotropic), <P is constant. A 
constant scalar field is the standard (simplest) way in the literature for causing 
exponential expansion, and in the isotropic case can be interpreted as creating 
a cosmological constant. For positive k, the solution approaches isotropy at 
late times. 

According to recent astrophysical data, the deceleration parameter of the 
universe lies in the interval —1.72 < q < —0.58 [3B]. In all models of this 
subsection at late time q converges to —I, so in this sense these models are 
consistent with the observational results. 



6.3 Shifted-reversed power seed 

We may also use u(t) = n ° 3 t ) m ; where CI3, £ c and m are non-zero constants, 
since this function is similar to the FRW Big Rip models. Then Eqs. (|4l?)) - (j5"l"j) 
become 

A(t) = -r-(t c - t) 1+2m [ha 3 3 (t c - t)' 3m + l 2 ] , (80) 
Bit) = — ^(£ c - £) 1+m / 2 [hal(t c - £)- 3 ™ + h] V2 , (81) 



a?m 2 (£ c -£)- 3m - 2 , . 

w Z 2 M(£c-£)- 3m + J 2 ] v ' 

The substitutions rn —> —n and £ t c — t bring these expressions into the 
same form as Eas. ([5T)) and Therefore all solutions from Section [BTTI apply, 
we just have to read the lines of the tables in that section from right to left, 
and switch the sign of n (and of course, shift any "i" designation to the other 
side of the singularity). So this seed does not only give Big Rip models (see 
Ref.[3"T] for detailed discussion). 



6.4 Other seeds and parameter choices 



Of course, v(t) being a free function, an infinite number of choices is possible, 
not to mention the freedom in the choices of l\ and li- Among these, we are 
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naturally interested in those choices that give simple metric functions A(t) 
and B(t). 

One way of achieving this would be to set h or I2 equal to zero. But I2 = 
takes us to isotropic solutions, which are not in line with the emphasis of the 
present work. The choice l\ = 0, although it will give simple-looking solutions, 
limits us to the case of massless field, and there seems to be no motivation 
for choosing a particular solution above others, except possibly v{t) — t (since 
the comoving volume stays constant), but this is already covered in subsection 

One can also choose combinations to make both v(t) and (his 3 + h) simple. 
For example, choosing v{t) = cos 2 / 3 (i), Iq = —2/3, l\ = —1 and I2 = 1 yields 

ds 2 = -dt 2 + cos 4/3 (t) [tan 2 (i) dx 2 + dy 2 + dz 2 ] (83) 

a metric with scale factors 

. , . sin(t) . 

A ^ - ^M) m 

B(t) = cos 2/3 (i) (85) 
a(t) = [sin(i)cos(f)] 1/3 (86) 

i.e. a set of universes that start with a 1-D Big Bang and end with a cigar-like 
Big Crunch (for ^ t ^ 7r/2, in appropriate units of time) or start with a 
cigar- like Big Bang and end with a 1-D Big Crunch (for n/2 ^ t ^ 7r). 
Another solution, 

ds 2 = -dt 2 + sinh 4/3 t (coth 2 t dx 2 + dy 2 + dz 2 ) (87) 

is mathematically similar, but describes universes with infinite lifetime, for 
example the t > universe starts with a cigar-like Big Bang and approaches 
de Sitter-like isotropically expanding universe at late time. 



7 Solutions in Einstein frame 

It is known that a scalar-tensor theory can be transformed to the so-called 
Einstein frame, where the gravitational scalar field becomes minimally cou- 
pled to curvature. The price to pay for this simplification is the equivalence 
principle: Massive point particles do not follow geodesies any more, in contrast 
to the Jordan frame, which we used so far in this work [HlfTB] . For the BD 
theory, the case W = —3/2 is special: We derived it from the vanishing of the 
Hessian (|25[) . it represents a fixed-point of the conformal transformation [15], 
and marks the boundary beyond which ghosts appear [37]. In the Einstein 
frame, this specialness is reflected in the gravitational scalar field becoming 
non-dynamical. Hence, it would be interesting to look at our solutions also in 
the Einstein frame. 
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The conformal transformation 

g ab = Q- 2 g ab , g ab 



2~ab 



takes us from the Jordan frame to the Einstein frame. Defining <P = e CT , the 
W = —3/2 BD action given in Eq. ((27)) takes the form 



A 



BV 



d x-s/—ge 



R + 2<J c a c 



EM*)) 



(89) 



where U\(o) = U ($)/<!>. It is easy to show that under a choice of a conformal 
factor 



the above action (1891) transforms into 



A £ = 



XyJ —g 



R-U 2 (<j)) 



(90) 



(91) 



This is exactly the Einstein-Hilbert action with a potential ^(c) = e a Ui(cr) = 
U(<P)/<1> 2 , and using the potential obtained in (|2"9")) it takes the form [^(ct) = A, 
where A is a constant and could be interpreted as the cosmological constant in 
Einstein frame. In this case, the (transformed) scalar field a does not appear 
in the action, hence the scalar field is non-dynamical as stated in the beginning 
of this section; and the field equations take the form 



Rab — -^9ab- 



(92) 



where the Ricci tensor R a b refers to the transformed metric g ab - 

The LRS Bianchi I spacetime can be brought back to its original form 



ds 2 = -dP + A 2 dx 2 + B 2 (dy 2 + dz 2 ) , 



(93) 



by a simple coordinate transformation after the conformal transformation. The 
transformations of time coordinate and scale factors from the Jordan frame to 
the Einstein frame are given by 



i = J V$dt, A = V$A, B = V&B. 



For this spacetime in Einstein frame the field equations (|92[) give 

2 



A B 
2-^ 

AB 

B" 
2-^- + 

B 



A 
2' 



A^_ 
~A 



►(f)- 

A'B' 
B AB 



A 
2" 



(94) 

(95) 

(96) 
(97) 
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where the prime represents derivative with respect to tilted time coordinate t. 
These equations can be solved exactly, giving 

A = c 3 (t + Cl )~ 1/3 , B = c 2 (t + Cl f /3 , (98) 

for A = 0; 

A = C3 sinh(fcf + ci) cosh _1 ^ 3 (fcf + ci), 

B = c 2 cosh 2/3 (fcf+ci), (99) 

for A = 8fc 2 /3 > 0; 

A — c 3 sin(fci + ci) cos~ 1/3 (fci + ci), 

B = c 2 cos 2/3 (fct + ci), (100) 

for A = — 8fc 2 /3 < 0. We note here that the obtained scale factors given in 
(|98p represent the well known Kasner solution. We would like to check the 
consistency between these solutions and those in the Jordan frame, (|4l?)) - ([5"Tj) . 
For given A, the Jordan frame solutions contain the arbitrary function v(t) 
and two arbitrary constants Iq and I2 (h is not independent). The Einstein 
frame solutions contain one arbitrary function <P(t) (which does not appear 
in the metric, however) and three arbitrary constants, c\, c 2 and C3. To show 
consistency, we need to transform the Jordan frame solutions to the Einstein 
frame and show that they agree with solution found in that frame. 

Applying the transformations in to (|4"9")) and ([5U|) . using (|5T|) . we get 



B = v{t) (102) 

Using (|102p in ((ST]) , we can write for <P in the transformed coordinates 

_ ^2 

but using l[9"4|). 

■~ dB dB dt ~. r- BB' 2 $ , . 

B = lit = S* = *^* = -<-*- ^ (104) 



11 (hB 3 +i 2 y 



so disappears from the equation: it cannot be determined in the transformed 
coordinates. This was to be expected, since the scalar field is absent from the 
action (|91|) . The arbitrariness (or information) in v(t) in the Jordan frame has 
shifted to $(t) in the Einstein frame. 

Returning to (|104[) and using the solution found for B{t) , e.g. for positive 

A, 

ll (ZiC2Cosh 2 (fc?+ci) +l 2 ) =c^-fc 2 sinh 2 (fcf+ci), (105) 
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where it should be noted that fc 2 = 3A/8 = 9Z x Zg/4 (see after ea.(fI5])). This 
will hold, if C2 is chosen such that I2 is equal to —he?,. This choice can be made, 
since obviously the set {ci, 02,03} cannot be independent of the set {/n,^}: 
If two solutions describe the same physical reality, the parameters (arbitrary 
constants) of one should be expressible in terms of the parameters of the other, 
although information could hide in the arbitrary function(s) in this case. 

Similarly, (|101[) will agree with if C3 is properly chosen; but c\ is not 
related to Zn or I2, it comes from the integration in the first term in f|94|) . 
The same calculations can be made for the cases of negative and vanishing 
A, establishing consistency of Einstein frame solutions with the Jordan frame 
solutions. 

For example, we can transform the solutions of subsection 6.1.1 (I2 = 0, 
isotropic) as 



3rj 
2kt 



A 



2kai 



B = ai t\e- 



(106) 



and those of subsection 6.1.2 {l\ = 0) as 

1/2 



n 2 a\ 



t (3/2)n- 



t = t Q (t + c x f'^ n \ A = A (t + c 1 )-^ 3 , B = B (t + c 1 ) 2 / 3 , (107) 



where 



tr 



4a? 



l/(3n) 



An 



2h 

3l 



1/3 



and Bn = 



1/3 



(108) 



Both solutions (|106p and (|107[l identically satisfy the Einstein frame equations 
(|95j) - (|97[) . if one recalls that l\ = implies vanishing of A. The constants t\ in 
(11061) and c\ in (|107[) correspond to the constants c\ that appear in ((98|) - (|100p . 

To summarize, the solutions we found in the Einstein frame, although con- 
taining only three parameters, are the transformed versions of the solutions in 
the Jordan frame, which contained an arbitrary function. This correspondence 
between a solution for the metric in the Einstein frame and an infinite number 
of metrics in the Jordan frame is possible, since the scalar field, which deter- 
mines the transformation between the two frames, is arbitrary in the Einstein 
frame. 



8 Concluding remarks 

In this paper we have examined the scalar-tensor Brans-Dicke theory of gravity 
for LRS Bianchi type I spacetime admitting Nocthcr symmetry. This symmetry 
approach is important because it provides us with a theoretical motivation 
to select a region of the solution space (For other motivations, leading to 
different regions of the solution space for the Bianchi type I spacetimes, see 
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e.g. Refs.|3H]-I12)- The Lagrangian density (fTij|) of LRS Bianchi type I becomes 
degenerate for W{<&) = —3/2, when the Brans-Dickc coupling function F(<P) = 
4> is used. This degeneracy is required for nontrivial solutions, hence we use this 
value as the BD parameter. The existence of Noether symmetry also restricts 
the form of potential U(<P), and allows us to find a transformation given by 
(jl2")l in which the metric potentials and the scalar field are stated in terms 
of new dynamical variables (/i, v, u), where the variable fj. is cyclic. Under the 
transformation (|42|) the Lagrangian (fT5)l reduces to a new, simpler one (|44|) . 

We have obtained the new set of field equations p5j) - (jT7)) for the LRS 
Bianchi type I spacetime by using these transformations. We have found the 
general class of solutions of BD field equations with potential U(<P) — \<P 2 in 
the background of LRS Bianchi type I spacetime exhibiting Noether symmetry. 
This solution family contains an arbitrary function, called v(t) in this work, 
and two arbitrary constants. 

In Section^ we gave examples using some simple forms of the seed function 
v(t) \ first as powers, then exponentials, powers of (t c —t), and some others. The 
first three arc chosen to be similar to popular models in FRW cosmology. The 
solutions are shown concisely in tables, clearly showing the relation between 
the behavior of the model universes and the parameters of the seed function. 

Because of the specialness of the value W — 3/2 for the BD-theory when 
one considers conformal transformations, we considered the problem also in 
the Einstein frame in Section [7] For this value, the scalar field becomes non- 
dynamical, taking on the arbitrary nature of the function v(t) in the Jordan 
frame. The solutions found in the Einstein frame are consistent with those 
found in the Jordan frame; they can be transformed into each other. 

The models found in Section [5] show a wide range of behaviors, featuring 
Big Bangs, Big Crunches, Big Rips, Bounces, various singularities of the cigar 
or pancake types, etc. While some of these solutions are of theoretical interest 
only, there are many expanding-universe solutions with acceleration, consistent 
with observational data. For example, all the solutions in Table [3] and all 
solutions in subsection l6.2l cxccpt the few strictly contracting solutions, feature 
late-time acceleration. 

Even the isotropic special case with power-seed function displays surpris- 
ing richness, despite the time-dependence of the scalar field being the same 
for all powers. One solution is particularly interesting: It is Minkowski space, 
containing a nontrivial scalar field. The particular nonminimal coupling and 
the potential selected by the Hessian Determinant condition and the Noether 
symmetry make this possible. Given the important guiding role of the concept 
of symmetry in modern theoretical physics, we believe that the family of solu- 
tions found and analyzed in this work constitute a potentially more relevant 
set among all possible solutions for LRS Bianchi typc-I cosmological models 
containing a Brans-Dickc field. 
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